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Deconfined regions in relativistic heavy ion collisions are limited to small volumes surrounded by 
a confined exterior. Using the geometry of a double layered torus, we keep an outside temperature 
slightly lower than the inside temperature, so that both regions are in the SU(3) scaling region. 
Deconfined volume sizes are chosen to be in a range typical for such volumes created at the 
BNL RHIC. Even with small temperature differences a dependence of the (pseudo) deconfining 
temperature on a colder surrounding temperature is clearly visible. For temporal lattice sizes 
N r = 4, 6 and 8 we find consistency with SU(3) scaling behavior for the measured transition 
temperature signals. 
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1. Introduction 

In relativistic heavy ion collisions (RHIC) we are dealing with an ensemble of volumes with 
typical length scales of 5 to 10 fermi. Particular volumes depend on details of the collision. These 
length scales are not really small when compared with a characteristic correlation length set by the 
inverse (infinite volume) deconfining temperature T c . Estimates of T c vary by up to 20% [1], but for 
our SU(3) investigation an approximate value is entirely sufficient, which we choose as in earlier 
work [2] to be 

r f = 174MeV (1.1) 

and hence T c 7 l = 1.1 3 fermi, 10% to 20% of a typical volume extension. Therefore, one can expect 
finite volume corrections due to cold boundaries, while with some notable exceptions [2, 3] most 
work in the literature addressed the infinite volume limit. 

In lattice simulations of pure SU(3) gauge theory equilibrium configurations at temperature 

T = WW, (12) 

can be generated by Markov chain Monte Carlo (MCMC) on 4D hypercubic lattices of size A T A S 3 , 
N t < N s with lattice spacing a and f3 g = 6/g 2 , where g is the bare SU(3) coupling. The infinite vol- 
ume limit is quickly approached through use of periodic boundary conditions (PBCs). Obviously, 
they are not suitable for our purposes and need to be replaced by boundary conditions (BCs) which 
model the confined phase. 

In [2] zero outside temperature defined by the strong coupling limit a(f5 g ) — > oo for j8 g — > 0, 
called cold boundary condition (CBC), was targeted. For realistically sized volumes corrections in 
the range from 30 MeV down to 10 MeV were found compared to signals of the deconfinement 
transition relevant for large volumes. Although agreement with SU(3) scaling was seen when 
varying N T from 4 to 6 and 8, one may be worried because the construction mixes a SU(3) scaling 
region with BCs from the extreme strong coupling limit. 

To use BCs from SU(3) configurations, which are also in the SU(3) scaling regions, the con- 
struction of a double-layered torus (DLT) was explored in [4, 5]. For the DLT boundaries are glued 
together as indicated for the 2D case by the arrows in Fig. 1 (the 2D topology is that of a Klein 
bottle). In the next section we discuss preliminaries for our investigation and in section 2 we report 
simulations using DLT BCs in 3D and PBCs in N t direction. We keep one DLT layer, called "out- 
side", below T c at T out « 160 MeV and search for signals of the deconfining transition in the other 
layer by varying its temperature. 



2. Preliminaries 

We use the single plaquette Wilson action on a 4D hypercubic lattice. Numerical evidence [6] 
supports that SU(3) lattice gauge theory exhibits a weakly first-order deconfining phase transition at 
coupling constant values j3/(A T ) = 6/g 2 c (N x ). The scaling behavior of the deconfining temperature 
is 

T c = c T A L (2.1) 
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Figure 1: Double layered torus in two dimensions. 



where the lambda lattice scale 

aA L = fl (pi) = Afe 2 ) (b g 2 )- bl/{2 ® e-"^ , (2.2) 

has been determined in the literature. The coefficients bo and b\ are perturbatively obtained from 
the renormalization group equation, 

11 3 J , 34 / 3 \ 2 

Using results from [7], higher perturbative and non-perturbative corrections are parametrized [8] 
by 

X(g 2 ) = l +e ln "'e- a ^ 2 +a 3 g 2 + a 4 g 4 (2.4) 

with In a i = 18.08596, a 2 = 19.48099, a 3 = -0.03772473, a 4 = 0.5089052. In the region accessi- 
ble by MCMC simulations this parametrization is perfectly consistent with an independent earlier 
one [9] (less than 1% deviation in the range of validity of [9]) and has the advantage to reduce for 
g 2 — > to the perturbative limit. 

For PBC finite size corrections to the deconfining temperature T c are negligible, so that we can 
use the N t = 4, 6, 8 estimates of the transition coupling j3|(A^) from [7] to fix Al- Using (1.2) 

Al = cT l f x (PI) = oT y f x (pf) aN,T c =f x m N T 174MeV = 5.07MeV (2.5) 

independently of the N t value. 

For the 3D DLT we use different coupling constants in the two layers. Along the lines of the 
usual derivations [2] the temperature in the entire DLT is given by T = \/[a(fii)N z ], x = in, out 
(1.2). We just have a system for which the coupling depends on the space position as it does in 
other systems like, for instance, spin glasses. The physical temperature does vary because the 
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I out 




158.15 MeV 




161.80 MeV 


N T 


4 


6 


8 


4 


6 


8 


Pout 


5.6500 


5.84318 


6.00458 


5.6600 


5.85514 


6.01821 


N s 


20 


30 


40 


20 


30 


40 


N s 


24 


36 


48 


24 


36 


48 


N s 


28 


42 


56 


28 


42 


56 


N s 


32 


48 


64 


32 


48 


64 


N s 


36 


54 




36 


54 




N s 


40 


60 




40 


60 





Table 1: Overview of our MCMC runs. 



lattice spacing depends on the coupling fi 8 = 6/g 2 . When using two different /3 g values, the lower 
one gives a lower physical temperature. We choose it to correspond to a temperature in the confined 
phase and call its layer "outside". The fi 8 value in the other layer, called "inside", will be iterated 
to find the maximum of the Polyakov loop susceptibility for its layer, which is our signal for the 
deconfining transition. To each plaquette we assign one of these values, fi 8 n and fi 8 ut , in a slightly 
asymmetrical way: If any link of plaquette is from the inside layer, fi 8 n is used, otherwise [} 8 uV So 
the inside lattice becomes slightly larger than the outside one. 

3. Numerical results 

In our simulations we accommodate both layers in the SU(3) scaling region, so that (2.2) 
holds, which starts slightly below j8 g = 5.65. With our smallest N t = 4 temporal extension this 
sets the lower limit on f5 8 ut . To work with round numbers we decided to start 7V T = 4 simulations at 
Afut = 5-65 and j3,f ut = 5.66. Translated into temperatures this corresponds to T out = 158. 15 MeV 
and Tout = 16 1.80 MeV. For larger N T = 6 and N T = 8 lattice these numbers have to be translated 
back into coupling constant values, listed in table 1 where we give an overview of the performed 
runs. The accuracy of the given numerical precision of the T out values reflects the SU(3) scaling 
relation when translating back and forth, but has no significance for the accuracy of the physical 
temperature. 

The small temperature difference between the confined phase and the infinite volume transition 
temperature makes it rather CPU time demanding to locate our deconfinment signals, the maxima 
of the Polyakov loop susceptibility. This is in essence overcome by brute force computer time. 
Relying on Fortran MPI code which is documented in [4], the bulk of our production was run 
at NERSC on their Cray XT4 and (more recently) Cray XE6, while code development and very 
limited production was done at FSU. All our runs used a statistics of 2 13 sweeps for equilibration 
and 2 21 sweeps for data production to overcome the statistical noise. Improved measurements 
techniques were tried, but did not result in major CPU time gains. Production on a 4 x 32 3 DLT 
used 1,024 processors and took 5.7 his in real time. Approximate run times on other lattices can be 
scaled from this. 
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Figure 2: Pseudo transition couplings from maxima of the Polyakov loop susceptibility. 



Fig. 2 plots the raw data for pseudo transition couplings J3| as determined from maxima of the 
Polyakov loop susceptibility and plotted versus the lattice ratio N T /N S together with polynomial fits 
in N r /N s up to order (N v /N s ) . For N T = 4 and N r = 6 we have excluded data from our smallest 
lattices, N T /N S < 0.5, because they suffer from peculiar finite size effects, which can be traced 
to the edges of the DLT layers. From down to up the lattices sizes are given by N T = 4, 6, 8. 
The corresponding j8 § ranges of pseudo transition temperatures shown in the figure increase from 
[5.69,5.725] (N T = 4) to [5.89,5.93] (N r = 6) to [6.06,6.11] (N T = 8). As seen in Fig. 3 the curves 
from these disconnected regions collapse into almost one curve when the SU(3) scaling relation 
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Figure 3: Pseudo transition temperatures versus volume edge length. 
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Figure 4: Pseudo transition temperatures versus volume edge length for CBC, DLT and PBC 
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(2.2) is applied to convert to T c in units of MeV, which we plot versus the volume edge length 
in fermi. 

Our final Fig. 4 presents T c results for outside temperature T out = 158.15MeV together with 
those of CBC from [2], where the strong coupling limit p 8 — > was used as BC. Even with our 
present high outside temperature the correction to the infinite volume critical temperature turns out 
to be remarkably large. So, this effect warrants further investigation. Instead of using a DLT one 
may incorporate cold boundaries by a surrounding region and, in particular, one may be interested 
in the inclusion of quarks in such calculations. 

4. Summary and conclusions 

• Even with a large outside temperature of about 160 MeV finite size correction of the SU(3) 
deconfining temperature due to small volumes are clearly visible. 

• Our estimates of (pseudo) transition temperatures show SU(3) scaling behavior when in- 
creasing the temperature extension of the lattice from N? = 4 to N T = 6 and N? = 8. 

• Small volumes increase the SU(3) transition temperature, while quark effects decrease it. If 
the increase holds also with quarks included, previous lattice calculations would underesti- 
mate the transition temperatures at RHIC. 
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